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Abstract 

The in-in effective action formalism is used to derive the semiclassical correction 
to Einstein's equations due to a massless scalar quantum field conformally coupled to 
small gravitational perturbations in spatially flat cosmological models. The vacuum 
expectation value of the stress tensor of the quantum field is directly derived from the 
renormalized in-in effective action. The usual in-out effective action is also discussed 
and it is used to compute the probability of particle creation. As one application, 
the stress tensor of a scalar field around a static cosmic string is derived and the 
backreaction effect on the gravitational field of the string is discussed. 







1 Introduction 



Our picture of the evolution of the early universe relies in the so called semiclassical theory 
of gravity which describes the interaction of quantum fields with the classical gravitational 
field. Order of magnitude arguments and Heisenberg's uncertainty principle tell us that 
there must be a period in the universe evolution, well after the Planck time, when the 
quantization of the gravitational field may be ignored but still the scale of its time variations 
is short enough to create elementary particles, so that matter quantization cannot be 
ignored. Since we lack a theory of quantum gravity it is still not known to what extent 
and in what sense this theory may be considered as a true semiclassical limit of quantum 
gravity interacting with matter fields. Plausibility arguments have been advanced by Hartle 
and Horowitz |IJ who show that the quantum corrections to the classical action of gravity 
interacting with N identical non self-interacting matter fields reduce in the leading-order 
1/N approximation to such semiclassical theory. 

The semiclassical approach provides the framework for some realistic scenarios which 
may explain some of the features of the present universe. One of these scenarios is in- 
flation which may explain the homogeneity and flatness problems of the standard big 
bang cosmology. In the inflationary model the quantum fluctuations of the inflaton field 
may be the source of the small gravitational inhomogeneities which seed galaxies or grav- 
itational waves. This may explain the universe large scale structure and the presence 
of a hypothetical background of gravitational radiation |4]]. Another scenario is the possi- 
ble formation of topological defects || as the universe undergoes some phase transitions. 
Topological defects, in particular cosmic strings, may seed structure || and may be an 
alternative to inflation for the generation of structure in the universe. 

In both scenarios the picture of the gravitational field that emerges is that of a confor- 
mally flat Friedmann-Robertson- Walker (FRW) background in which small gravitational 
perturbations are present. Large anisotropies and inhomogeneities might be present only 
if the universe had emerged highly inhomogeneous from the Planck era into the classical 
regime J7[. But since the quantum consequences of a highly inhomogeneous model are diffi- 
cult to estimate one assumes that by studying small perturbations on a FRW background, 
a qualitative picture of the evolution of the more extreme case may result. 

Here we are interested in the quantum effects produced by the presence of small per- 
turbations in conformally flat backgrounds. Quantum effects due to small anisotropies 
were first considered by Zeldovich ||, Zeldovich and Starobinsky ||, Hu and Parker |I0[ , 
Hartle and Hu ||11|| , and Birrell and Davies who computed the creation of conformally 
coupled particles interacting with the anisotropies. Conformally coupled particles are not 
created in conformally flat backgrounds (FRW) ]13[ but the anisotropies break the con- 
formal symmetry. Different techniques were used for such computations, these techniques 
go from a perturbative evaluation of the Bogoliubov transformations relating two vacua 
of the quantum field, to the evaluation of the in-out effective action of this field in the 
given gravitational background. These results were extended to the presence of arbitrary 
perturbations, including inhomogeneities, MM by a technique based in the perturbative 
evaluation of the scattering matrix which had been used in flat backgrounds by Sexl and 
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Urbantke [|HJ and Zeldovich and Starobinsky |I6 

Quantum effects on the geometry, the so called backreaction effect, are more difficult 
to evaluate because this requires, on the one hand, the computation of the renormalized 
stress tensor of the quantum field in order to modify the classical Einstein equations, and on 
the other hand, it requires the solution to these semiclassical equations. It was argued by 
Zeldovich || that the backreaction would tend to dissipate the inhomogeneities as in a sort 
of gravitational Lenz's law effect. This is a mechanism to homogeneize the universe, but it 
is usually not advocated because in the standard scenario one assumes cosmological models 
which cannot explain the present large scale homogeneity by any causal mechanism after 
the Planck era; the inflationary scenario, on the other hand, seems to solve the homogeneity 
problem quite naturally. It is nevertheless a mechanism for entropy production. Early work 
on the backreaction effect on the geometry due to anisotropies was done by Lukash and 
Starobinsky ]17 and Lukash et al. [18], who assumed very special conditions near the 

who considered a Bianchi type I anisotropic 
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Planck time, and by Hu and Parker 
model, evaluated the stress tensor in the low frequency approximation and computed the 
resulting modified Einstein's equations numerically. The results of such work indicate that 
the dynamical mechanism of particle production achieves a rapid damping of the anisotropy 
if the calculations are extrapolated to the Planck era. 

The computation of the quantum stress tensor is generally difficult in practice. How- 
ever, for small perturbations on a conformally flat background one may use perturbative 
methods to get explicit expressions. One of the most powerful and efficient methods and, 
one that is very well adapted to a perturbative scheme, is based on the one loop order 
computation of the so called in-in effective action for quantum fields interating with the 
gravitational perturbations. This technique is an effective action technique adapted to 
compute expectation values of quantum operators. It was first proposed by Schwinger JTU 
and Keldysh |20| and developed by Chou et al. plfl . Jordan [^2| and Calzetta and Hu 
p3j developed the technique on a curved background, and it was then applied to derive 



the stress tensor of a quantum scalar field coupled to small anisotropies on a cosmological 
background f23|. 

The use of effective action methods in the backreaction context was first considered 
and Fischetti et al. 



by Hartle [24 
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111 studied the effect of 



and Hartle and Hu 

anisotropies. But in their formalism the basic element is the usual effective action which 
is related to the generating functional of the in-out vacuum persistence amplitude. This 
in-out formalism leads to matrix elements rather than expectation values for the quatum 
operators. Thus, One does not get directly from the in-out effective action the vacuum 
expectation value of the stress tensor of the quantum field, and one still needs to compute 
the Bogoliubov transformation between the in and out vacua. This method is, however, 
very useful for the computation of the particles created, since the probability amplitude 
for particle creation is directly related to the vacuum persistence amplitude. 

In this paper we compute the in-in effective action to the one loop order for a massless 
scalar field conformally coupled to small gravitational perturbations on a spatially flat FRW 
background. The in-in effective action is used to derive the quantum stress tensor and the 
corresponding semiclassical Einstein's equations. Our results generalize the Calzetta and 
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Hu |23| results to the case of arbitrary small perturbations including inhomogeneities, and 
the stress tensor we derive coincides with that obtained by Horowitz and Wald [2fj], who 
used an axiomatic approach to derive it, and by Starobinsky [27] who used a modified 
Pauli-Villars regularization method @. One should stress, however, that the stress tensor 
computed does not include the energy of the particles created which is a second order 
correction to the computed terms; it includes only first order vacuum polarization effects. 
One might wonder that although the energy of the particles created is small it might have 
a long term cumulative effect. 

We should mention that the axiomatic approach to derive the stress tensor has been 
quite succesful in several situations. Thus Horowitz p8| obtained the stress tensor due 
to a scalar field minimally coupled to arbitrary gravitational linear perturbations on a 
flat spacetime background applying the axiomatic arguments outlined by Wald [p9| . This 
tensor was rederived by Jordan |HJ using the in-in effective action method; note that the 
case of conformally coupled fields may also be obtained from our cosmological model when 
the conformal factor is taken constant. Another approach to the quantum stress tensor 
based on an iteratively evaluated mode decomposition was developed by Davies and Unruh 

In this paper we shall not consider the solutions to the semiclassical equations, except 
in a simple example involving a cosmic string. The correct approach to this problem is still 
controversial. In fact, the semiclassical equations are known to admit runaway solutions as 
a consequence of the fact that they are dynamical equations with higher order derivatives. 
Horowitz [^| and Jordan found from these solutions that flat space is unstable against 



quantum effects. Whether these solutions are physical and thus signal a true instability, 
or unphysical and thus spurious, has been the subject of some discussion in recent years. 
Simon [32| has argued that the semiclassical correction to Einstein's equations must be 
seen as analytic perturbations, in terms of the Planck constant h, to the classical Einstein's 
equations and that, as such, only solutions which are also analytic in h are physical. A 
consistent perturbative approach to find reduced equations, i.e. dynamical equations which 
are second order at each order of perturbation is known |33fl . When this is applied it is 
found that flat space is perturbatively stable to first order in h ||32|| . Reduced semiclassical 
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equations have been obtained also in some cosmologies [34 



35 1, on the other hand, 



has argued, on the basis of how the stress tensor is renormalized, that this tensor cannot 
be considered the first term of an expansion and therefore the previous reduction methods 
should not be applied. 

In order to make this paper reasonably self-contained the in-in effective action formula- 
tion is summarized in section 2 with a view to practical applications. In section 3 the in-in 
effective action to the one loop order is derived for a scalar field conformally coupled to a 
nearly conformally flat metric. Since along this computation one derives also all the terms 
needed for the in-out effective action, this action is also discussed, and it is used to derive 
the probability for pair creation; the results agree with those obtained by other methods. 
In section 4 the stress tensor for the quantum field is derived from the in-in effective action 
and the semiclassical correction to Einstein's equations is written down. The stress tensor 
is seen to agree with that obtained by Horowitz and Wald [|2B| and Starobinsky |27 



As 
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an exercise the semiclassical equations in two dimensional spacetime are also derived using 
the same formalism. In section 5 we apply the previous formula to compute the stress 
tensor of a quantum field around a static cosmic string and we discuss the backreaction 
effect on the gravitational field of the string. Note that since the gravitational field of 
a cosmic string can be considered a small perturbation on a flat background the above 
perturbative technique (in the sense of metric perturbations) can be applied. The results 
are in agreement with those found by other non perturbative methods f36fl ||37|| but this 
perturbative method opens the possibility of computing the quantum stress tensor even in 
time dependent situations. Work along these lines is in progress. 



2 In-in functional formalism 

In this section we summarize the in-in functional formalism for the evaluation of the in-in 
effective action with a view to the applications of this paper. We follow, essentially, the 
presentations by Jordan fl2"2 |, Calzetta and Hu f23f , and Paz ]38 . 



Quantum corrections to a classical field theory can be studied with the help of the 
effective action. For simplicity, we consider the quantization of a scalar field <fi{x). The 
usual in-out effective action is based in the generating functional W[J] which is related to 
the vacuum persistence amplitude in the presence of some classical source J(x) by 



iW[J] = (0,out\tn,0)j. (2.1) 



This functional carries all the quantum information of the connected graphs of the theory. 

When one couples an external field J{x) it is convenient to use the interaction picture in 
which the states \ip) evolve in time according to the Schrodinger equation Hi\ip) = idt\ip), 
where Hj is the interaction Hamiltonian operator Hj = J d n ~ l xJ(x)(j)(x), <fi{x) is now the 
field operator in the Heisenberg representation and n the number of spacetime dimensions. 
The solution of this equation may be formally written as, 



Ita 



T®exp[i I' 2 dtHj) \^) tl , (2.2) 



where is the usual time ordering operator, and (2.1) can be written as 

jw[j] = ^ out \ T (t) exp (a J°° dtH ^ | i7lj o). (2.3) 



e 



It is easy to see from the classical field equations for <f>(x) in the presence of J(x) that 
expzW^fJ] satisfies the integro differential Schwinger-Dyson equation, and that one may 
give a path integral represenation for its solution as, 

e iW[j] = J p^ e i(SM+j« j (2 _ 4) 

where S[4>] is the classical action of the field theory and the common shorthand notation 
J0 for the integral / d n xJ(x)(f)(x) has been used. The functional integral is taken with 
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the following boundary conditions: <fi — > e T%ujt , where ui > 0, when the time t — > ±00, i.e. 
the scalar field has only negative and positive frequency modes in the in and out regions 
respectively; the interaction is assumed to be switched off at these asymptotic regions. By 
differentiating with respect to the source one generates matrix elements from W[J] , 

5W[J] _ (0,(mt\<Kx)\m,0)j _- 

SJ(x) {0,aut\in,0)j ~ n J ' 1 ' 

If we assume that the above expression can be reversed, the effective action is defined as 
the Legendre transformation of the generating functional, 

= W[J] - J0. (2.6) 

This functional of is the generator of the one-particle-irreducible graphs (graphs that 
remain connected when any internal line is cut) and contains all the quantum corrections 
to the classical action. From (2.6) one may derive the dynamical equation for the effective 
mean field 0[O], i.e. the matrix element of the field in the absence of the source J(x), as 



= 0, (2.7) 

=0[O] 

which expresses the quantum corrections to the classical equation as a variational problem 
of the effective action. 

In order to work with expectation values rather than matrix elements one can define a 
new generating functional whose dynamics is determined by two different external classical 
sources J+ and J_, by letting the in vacuum evolve independentenly under these sources, 

e iW[j + ,j.] = T£(0,in\a,T)j_(a,T\in,0) J+ . (2.8) 

a 

Here we have assumed that {\a, T)} is a complete basis of eigenstates of the field operator 
4>{x) at some future time T, i.e. x)|a, T) = a(x)|a,T). Then (2.8) may be written 
according to (2.2) as 

e iW[J + ,J-] = J da^inlT^e-^-^^^-^^l^T) 

x(a,T\T^e l ^ dt ^ dn ' lxJ+{xmx) \ in ,0)^ (2.9) 

where T® and mean, respectively, time and anti-time ordered operators and da means 
da = n x c?a;(x) where x are the points of the hypersurface £ defined by t = T. The 
generating functional has also a path integral representation, 

e iW [J+ ,j.] = JdaJ £>[</>_] e -^[4>-]+J-4>-) J P[ + ] e i(5[0 + ]+j + + ) ) (2 1Q) 

with the boundary conditions that + = 0_ = a on X and that the fields + and <p_ are 
pure negative and pure positive modes, respectively, in the in region i.e. <p± — > e ±lult at 
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t — > — oo (vacuum boundary conditions in the remote past). In a more compact form one 
may write, 



JW[J+,J- 



V\ 



i{S[<p + ]+J + d>+-si4>-]-j^4>-} 



(2.11) 



where it is understood that the sum is over all fields </>+, 0_ with negative and positive 
frequency modes, respectively, in the remote past but which coincide at time t = T. These 
boundary conditions can be made explicit by substituting m 2 by m 2 — ie, where m is the 
field mass, in S[<p + ) and by substituting m 2 by m 2 + ie in S [</>_]; the latter is also some 
times indicated by writing S*[<f)J\ instead of S[(j)-] [p3| . This integral can be thought of as 
the path sum of two different fields evolving in two different time branches [ 3~9fl , one going 



forward in time in the presence of J + from the in vacuum to a time t = T, and the other 
backward in time in the presence of J_ from the time t = T to the in vacuum, with the 
constraint <p + = 0_ on E. Because of such a path integral representation, this formalism 
is often called closed time path formalism. 

The functional W[J+, J-\ generates expectation values of the field rather than matrix 
elements. For instance, we have, 



6W[J+, J_ 



5J 4 



(O,in\<j)(x)\in,O)j = 0[J}, 



(2.12) 



j±=j 



instead of equation (2.5). This functional generates not only the desired expectation values 
of time ordered field operators but also the anti-time ordered ones in the same footing 

j: iW[J+,J.] 

(0, inlTMMjfc) • ■ ■)T^( ( j ) (x 1 ) ■ ■ -)|m, 0). (2.13) 



iSJ^xt) ■ ■ ■ (-i)5J-(y 1 ) 



j ± =j=o 



In analogy with the in-out formalism the in-in effective action is defined as the Legendre 
transform of the new generating functional as, 



r[0 + ,0_] = wqj + ,j_]- j + l 

where the external sources are functionals of the fields 

SW[J+, J_] 



5J 4 



±0 ± [J + ,J_ 



+ J-0-, (2.14) 
) + and 0_, through the definitions 

(2.15) 



which we assume can be reversed. 

From the definitions (2.14) and (2.15) we get the equation for the expectation values 
4>±[J+, J-]j, i.e. 



(2.16) 



and by taking J± = in (2.15) we recover the equation for the vacuum expectation value 
of the field 0[O] = 0±[O,O] = (0,in\<f>(x)\in, 0), 



0. 



(2.17) 



0±=^±[O,O]=0[O] 



6 



This equation does not follow from a simple variational principle in terms of a single field 0: 
in the in-in action we have two fields 0+ and 0_ which are treated independently and only 
when the sources have been eliminated they become the vacuum expectation value. Note 
also that r[0, 0] = as a consequence of (2.14) and of the fact that W[J, J] = 0, which 
follows from (2.8) and the usual normalization for the states. Equation (2.17) is a dynamical 
field equation which admits an initial value formulation: the solution 0[O] = <f)±[0, 0] is real 
and causal, i.e. the solution at one spacetime point depends only on data on the past of 



that point 22 



For a free field theory, i.e. a theory with a quadratic action, we can compute W / o[</+, J-\ 
from (2.11), which becomes now a Gaussian integration for the two independent fields 
+ and The corresponding propagators will be determined by the very particular 
boundary conditions of this problem. In fact, let us assume that the free action for + is 



Sole 
for 



/ c? n x|(9 M + 9 M + + (m 2 — ie)4> 2 + ) and that we have an analogous action Sq[ 



then the classical field equations are, 

(□ - m 2 ± ie) 



-J±(x) 



(2.18) 



At this stage we can introduce the compact notation, 



<t>a{x) 



S[<f> + ] - S* 



Jaix) 




(2.19) 



where a and b take the two values + and — , to simplify the mathematical expressions. The 
solutions of the classical field equations, that satisfy the boundary conditions, <p° ± — > e ±luJt , 
when t — ► — oo with uj > 0, and 0+(T, x) = <p°_(T, x) in the hypersurface E, which we take 
here at t = T — > oo, can be written as, 



<fa{x) 



d n yG° ab (x,y)J b (y), 



where G® b (x,y) is the matrix 



G° 



ab 




(2.20) 



(2.21; 



defined with the Feynman, A F , Dyson, A D , and the positive, A + , and negative, A , 
Wightman functions: 



A F (x - y) 
Az?(x - y) 



d n p e 



ip-(x-y) 



A ± (x - y) = (±2vri 



(2ir) n p 2 + m 2 — ie ' 

d n p e ip - {x ~ y) 
(27r)™ p 2 + m 2 + ie ' 

d n p 



(2.22) 



(27T) 



e ^-y)S(p 2 + m 2 )8{TP°). 
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These propagators are solutions of the equations 



A ac G° cb (x,y) = 5 n (x-y)5 ab , (2.23) 
where the operator A° ab is the diagonal matrix defined by 

A° ab = diag [(□ - m 2 + ie\ , - (□ - m 2 - iej] . (2.24) 
Furthermore the Feynman and Dyson Green functions have a mode decomposition 
A F (x -y) = 6{x° - y°)A-(x - y) - 6(y° - x°)A + (x - y), 

(2.25) 

A D {x -y) = 9(x° - y°)A+(x - y) - 9(y° - x°)A"(z - y), 

which reflect the boundary conditions imposed over each classical field solution <Pa{x), be- 
cause the Green functions A ± (x — y) ~ e ±iu-(x -y ) correS p 0nc l to negative and positive 
frequency modes respectively (here u = \/p 2 + m' 2 ). Note that it is also satisfied that 
dt(j>+(T, x) = d t (p°_ (T, x) at S. With these propagators to guarantee the boundary condi- 
tions the Gaussian integration of (2.11) for a free field is, 

W [J a ] = -i I d n xd n yJ a (x)G ab (x - y)J b (y), (2.26) 

where a term independent of J a has been discarded to satisfy Wq[J, J] = 0, and one can 
use now (2.13) to generate time ordered and anti-time ordered expectation values of field 
operators. In particular we have that 

(0,in\T®<f>(x)<f>(y)\in, 0) = iA F (x-y), 

(0,in\TW<f>(x)<f>(y)\in,Q) = -iA D (x-y), (2.27) 

(0, in\(j)(y)<p(x)\in, 0) = —iA + (y — x) — iA~(x — y). 
For interacting fields one can proceed as usual by writing 

e iW[j a ] = e ifd"xc mt (-^) e iWo[j a ]^ ^2_ 28 ) 

where we have separated the Lagrangian into a free and an interacting part, C = C Q + C int . 
Then one may continue in the usual perturbative fashion; however, we are not going 
to consider self-interacting theories in this paper, the only interaction will be with the 
gravitational field. 

Let us now proceed to the main objective, namely, the evaluation of the effective action 
T[0 a ] up to the one loop order, which corresponds to the first order expansion of VK[J ] 
in powers of h. As usual if we assume that the action is bounded from above then 
we can go to Euclidean space and solve (2.11) by the steepest descent method; we keep, 



8 



however, the Minkowskian notation. Let us denote by <f>+(x) and <f)^\x) the solutions of 
the classical field equations which may, or may not, include self-interactions, 



5S\ 



<.(°)i 



(2.29) 



and let us expand the exponent in (2.11) about these background fields: 
<S[0 Q ] + f d n xJ a (x)Mx)=S[<P { a ) } + / d n xJ a {x)^\x) 



+ - / d n xd n y </> a (x) 



X 



A ab (x,y) Mv)-<f>b\v) +"(2.30) 



where 



5cf> + (x)5(j) + (y) 



<t>+=<t> 



, A—(x,y) = - 



5 2 S* 1 



f2.3i; 



and, of course, A + _(x, y) = A- + (x,y) = 0. Substituting this into (2.11) the integration is 
now Gaussian and we can write to this one loop order, 



JW[Ja 



iw(0) ^ (detA ab (x,y))-^ 2 



(2.32) 



where W<® [J«] = S[4>^} + J d n xJ a <p^\ In terms of the propagator G, which is a functional 
of the background fields <pf\x) and takes a 2 x 2 matrix form , i.e. G(x,y) = A~ 1 (x ) y) 1 
we can write (2.32) as 



W[J a )~WW[J a ]-~Tr(\nG). 



(2.33) 



The effective action, which is a functional of <fi a , can now be explicitly found to the 
same order. Using (2.14), (2.15), and the fact that 4> a differs from <%°) by a term of order 
h we can show that W' '[J o ] — S[(j> a ] + / d n x.J a (j) a , so that finally we have 



r[0J ^ 



s\ 



-Tr (In G). 



(2.34) 



Now the equations for cf) a can be deduced from (2.16) using the explicit functional 
dependence on the fields given by (2.34). However we should note from (2.17) that in 
order to get the field equations for the expectation value of <f>(x) we only need the explicit 
dependence of the effective action on one of the fields, <f) + or </>_. Therfore we are only 
interested in the dependence of (2.34) on cf> + , say. Following Paz |38| we can write 



5<f) + (x) 



(Tr (In G)) 



dydzG ab (z,y)- 



6<p+(x) 
6 

dydzG ++ (z, y) - G++(z, y) 



6<j)+(x) 



5 

5(j)+(x] 



(Tr(\nG ++ )) 



(2.35) 
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where we have used that G ab = A a b is diagonal, see (2.31). Thus we have, 

~ S[j+] - l -Tr(\nG ++ ) + F (2.36) 



where F includes all the terms which do not contribute to the variation of the field <f> + , 
i.e. g^j^j - - — 0. This expression is very similar to the one loop in-out effective action 

which is given by the above equation (2.36) where C++ is substituted by the Feynman 
propagator Ap, the main difference is in the boundary conditions: the propagator G ++ is 
defined as an expectation value and not as an in-out matrix element. 

This formalism can be extended to curved spacetimes without difficulties assuming 



that the spacetime is globally hyperbolic [^2|. The hypersufaces of constant time are now 



Cauchy hypersurfaces and the in and out states are defined in the Cauchy hypersurfaces 
corresponding to the far past and far future respectively. Now the spacetime integrals must 
be performed with the volume element d n Xy/—g where is the spacetime metric. The 
above expressions (2.20), (2.23) and (2.36) are still valid except that now the Feynman, 
Dyson and Wightman functions have a different representation to that of (2.22). If the 
spacetime is asymptotically flat in the in and out regions the previous boundary conditions 
for the fields <f)±(x) will also apply; if not, in order to be able to define physically meaningful 
in and out vacua we must assume that we are still able to define positive and negative 
frequency solutions in the asymptotic regions. This is always possible, for instance, if the 
asymptotic regions admit approximate timelike Killing fields. But, generally, in a curved 



spacetime the in and out vacua are not equivalent. Jordan [22| has shown that for quantum 
scalar fields in a curved spacetime the field equations are real and causal up to the two loop 
order and he has also checked the unitarity of the formalism restricted to vacuum states. 

Before ending this section let us rewrite (2.36) in a more convenient form for us. In gen- 
eral the propagator G ++ cannot be found exactly and has to be evaluated perturbatively. 
For instance, in the next section we will take perturbations h^ u (x) to a given background 
metric and only the exact propagator corresponding to the background is known. Thus we 
write 

A ab = A° ab + (V^ + V^ + ...) (2.37) 
where A° ab is the unperturbed (diagonal) operator whose propagator, G° ah1 is known, 

KA = <U (2.38) 

and the diagonal operators V^ b + Vj^ + • ■ ■ contain the perturbative terms (in the next 
section they will correspond to perturbations of order \h^ u \ and \h^ v \ 2 respectively). We 
can write 

Gab = G° ab - G° ac (v}^ + ■ ■ -)G db 

_ n o rio T/(i)/nro r 10 \f( 2 )n® j_ r<o x/(i)r»o x/W/nro , / 9 o \ 

— ^ab ~ <^ac V cd ^db - ^ac V cd ^ db + (j ac V cd U de V ef U fb H 5 
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where the products are operator products. In particular, 

Lt ++ - °++ ~ Lr +a V ab Lr b+ ~ Lr +a V ab Lr b+ + Lr +a V ab ^bc V cd H J l^.4UJ 

expanding the logarithmic term in (2.36) and using that A° ++ G° + _ = 0, we finally get, 



W-] * S[0 + ]-^Tr(lnG° ++ )+F 



+^G%_F_ w G t i + + • • • J . (2.41) 

where we have defined V® = V+\ and = — ViL following (2.31). Note that if it 
were not for the last term which involves the propagator G° + _ this expression for the in-in 
effective action would agree with the in-out effective action which involves only one field 



4>(x); see, for instance, Hartle and Hu [11|. Therefore the term containing the propagator 
G+_ is the only new term that contributes to the field equation for 4>+(x). It can be seen 
38| that the effect of the last term in (2.41) is to make the field equation for 4>(x) causal: 



if one takes the derivative of the in-in effective action with respect to <f>+(x) and puts 
= 0__ = 0, the resulting field equation is causal. 

Notice that in the case of a free scalar field on a nearly flat background the propagator 
G° ab which solves (2.38), corresponding to the flat background, and which provides the 
boundary conditions of the in-in problem is simply given by G° ++ = Ap, G°__ = —Ad, 
G°_ + = A~ and G° + _ = — A + , as can be seen from (2.23). In fact, from (2.31) A is, in this 
case, the operator defined in (2.24). 



3 In-in effective action 

In this section we compute the in-in effective action (2.41) for a conformal field in a 
nearly conformally flat spacetime. The cosmological background consists of a spatially 
flat homogeneous FRW with small perturbations as, 

9ni>( x ) = a2 (v)(v^ + h^ u (x)) . (3.1) 

where a(rj) = exp(u(rj)) is the conformal factor, rj is the conformal time drj = dt/a, t 
is the cosmological time, and h^x) is a symmetric tensor representing arbitrary small 
perturbations; we take the metric signature ( — !-■••+). The classical action for a free [i.e. 
with no self- interactions) massless conformally coupled scalar field $(x) is given by 



S m [g^u, <l>\ = -\j d^ypti [r v d^d u <5> + i{n)M 2 ] , (3.2) 

where £(n) = 4 ^~ 2 ^ , R is the Ricci scalar for the metric g^, and we take the spacetime 
dimensions, n, arbitrary for the moment in view of dimensional regularization. Because of 
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the conformal coupling one may simplify the problem by defining a new matter field (f>(x) 
and a new metric g^ v (x) as, 



[x 



) = e^^^x), g^{x) = Vlu , + h^x). (3.3) 



Then the action (3.2), after integration by parts and assuming no contributions of the 
surface integrals, is equivalent to the action for the field (f>(x) in the nearly flat metric g^ u , 

S m [g,u, <t>\ = -\j d n xV=g [g^d^dA + £(n)it>0 2 ] , (3.4) 

where R is the Ricci scalar for the metric g^ v (x). Therefore the problem has been reduced 
to that of a scalar field <f>(x) in a nearly flat background. We must keep in mind that 
the physical field is $(x), but the fact that the two fields differ by just a power of the 
conformal factor, which is a function of time only, simplifies considerable the connection 
between the vacua of the two fields. For instance, a positive frequency mode in flat space 
will correspond to a positive frequency mode in the conformally related space. Since for a 
free field in flat space the in and out vacua are equivalent (there is no particle creation) the 
same is true for the vacua of the conformal field in the conformally flat background. As a 
consequence non trivial quantum effects can be produced only by the breaking of conformal 
flatness which in this case is due to the coupling of the quantum field with the gravitational 
perturbations. The above action can be expanded in terms of these perturbations as, 

oo 

S m [V» 0] = Si O) [0] + £ &h„ u , 0], (3.5) 

n=l 

where the first term is simply the action for the field <fi in flat spacetime, and the higher 
perturbative terms carry all the information on the interaction with the perturbations. 

Since we are interested in deriving the semiclassical correction to Einstein's equations 
due to the quantum effect of the scalar field but keeping the gravitational field as classical, 
we have to add to the effective action T m (m stands for matter fields), the classical action of 
the gravitational field S g [g^ u ]. We should also add the action of any other classical source 
but we shall ignore this for simplicity; note that its effects on the semiclassical equations 
may be taken into account by simply adding the corresponding classical stress tensor to 
the quantum stress tensor. Furthermore, in order to renormalize the effective action it is 
sufficient to add to the usual Einstein's action, terms quadratic in the Riemann tensor, 

s g [%A = /^(-^)) 1/2 {r^^) 

'R^{x)R^{x) - R^(x)R^(x)] } , (3.6) 



/i"" 4 

+ 



28807r 2 (n - 4) 

where \i is an arbitrary mass scale which will be useful in dimensional regularization. The 
quadratic terms with poles at n = 4 are those which are necessary to cancel the divergencies 
of T m ; notice that with this election one obtains the correct trace anomaly. 
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Following the previous section we introduce now two fields <f>+(x) and 4>-{x) which at 
some future hypersurface X coincide: <$+{x) = Since the scalar field is proved by 

the gravitational field we must assume that the two fields evolve in two different geome- 
tries g^ u (x) and g^ u (x) respectively where g^ v (x) = rj^ + h^ u (x), and since we assume 
that the fields have no interaction other than the gravitational we need not introduce the 
classical external currents J±(x). Thus we shall write the total in- in effective action for 
the gravitational and the matter fields as, 



r (ii)k h %\ = S 9 k Ku\ - S g [u, h ] + T m [h + 0+; h <p_] 



(3.7) 



where T m [h^ u , + ; h , </>_] contains the quantum effects of the scalar field. 

Following eq. (2.41) we write T m in a perturbative expansion in h^ v as T m = F$ + 
T$ + r£P + ... and write only the terms which contribute to the variation of h+ u , which are 



also those which contribute to the variation of <fi + (x). Thus we can write, V 



r 



(i) 
(it) 



SV> + T« and r 



Sg = 5} 2) +r« where 
= *W 



(o) 

(«) 



~Tr(lnG° ++ ), 



+ l -Tr{vi 2) Gl + ) 



Tr{v£ ) (? ++ V£ ) GP ++ ) + Jrr(y| 1) G°_yi 1) G' _ + ). 



(3.8) 



To compute V {1) and \/ (2) we use (2.31) and (2.37) and expand S. 



as in (3.5), 



h^d^d u <p + f (n) (V 2) + 



where /i 



2 hh^ii/ 



(3.9) 

(3.10) 
(3.11) 

2) 



X /i 2 ^ - iKph^r]^ and and R ( 
are the first and second order terms, respectively, of the scalar curvature (see equation B.9 
from the appendix). From these expressions one gets the operators A , and by 
diffenciation with respect to the field <f>(x), 



A°(f)(x) = n<f>(x), 

VW(x)<j>(x) = - [(d^ix)) d v + h^{x)d lx d v + an)R {1) ] <P(x), 



V {2 \x)<p{x) = {dji» v {xj) d u + h^(x)d^d u - £(n) [R {2) + -hR {1) 



<t>{x) 



(3.12) 
(3.13) 
(3.14) 
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The explicit form of the operator will not be needed, however. Now we can write the 
propagator, from (3.12), as, 

r d n k p ik -( x ~ x ') 

Gl + (x, X ') = A F (x,x') = -j (27r)ra k2 _ %e , (3.15) 

recall that this is the propagator for the field 4>(x), the propagator for the physical field 
$(x) is related to this by the conformal factor |L^], but we do not need it here. All we 
need to know is that the boundary conditions for the physical field are determined by the 
boundary conditions of <fi{x) in flat space. 

The effective action depends on the fields <j>± and the metric perturbations h^ u . We may 
obtain the equation for the field <f>(x) (field equation) by functional derivation with respect 
to + , i.e. equation (2.17), and the equations for the metric perturbations (backreaction 
equations) by functional derivation with respect to ht v . Our primary interest is to obtain 
the backreaction equations. From equation (2.17) one can see that the vacuum expectation 
value = (0, in\4>\in, 0) = 0{h^ u ) (it would be zero in flat space), and thus the contribution 
to the backreaction equations coming from the matter action term will be at least of second 
order in the perturbations when the field equations are substituted. Formally, one can 
compute the effective action as a functional of h^ u only, and thus, the terms involving the 
field (p + (S!$[(j) + ]) are not necessary. 

We can now compute each of the (divergent) terms (3.8). As it is well known the 
first term T$ is easily renormalized by adding a suitable counterterm that cancels the 
divergencies which lead to the conformal anomaly [12|] f25f , but, this term does not play 



any role in the stress tensor of the field because it is independent of h^ u . 

We can now go to the next term in (3.8), r$, but this formally divergent (tadpole) 
term has no contribution, since Tr(V^ G° ++ ) involves n-dimensional integrals of the form 
1/k 2 , k a /k 2 , and k a kp/k 2 (where k a is the integration momentum variable) which are 



identically zero in dimensional regularization | 4T| . Therefore there is no term linear in 
in the effective action, i.e. we have T$ = 0. 

The first non trivial quantum contributions to the stress tensor coming from T m are 
quadratic in h^ u and we thus need to compute T^. Here again the second (tadpole) term 
in gives no contribution in dimensional regularization, since the typical 

integrals are of the same type that those of Tr(V+^G+ + ). For this reason we do not need 
the explicit form of the operator V^. 

The problem is thus reduced to the evaluation of the third and fourth terms in (3.8). 
As we have stressed in section 2 the third term also appears in the evaluation of the in-out 
effective action and the fourth term is typical of the in-in contribution to this order. Let us 
evaluate the third term, i.e. T\ = -{Tr{V^ ] G° + + V^ ] G\ + ) (recall that G° ++ = A F , with 
m = 0), 



Ti = ~y d n xd n x'vi 1) {x)A F {x,x')v} 1) {x')A F {x',: 



4 J J (27r) n (27r 
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X 



iq-(x~x') 



X 



((d a ,hf{x')) d p , + hf(x')d a ,d^ + £(n)i4V 



e ip-(x'-x) 

p 2 — ie 



. (3.16) 



We now introduce the projector P^ = rf v — p^p u /p 2 , the symbol ?^ a/3 ', h p,v = i]^ aul3 h al 3, 
change the p integration by p' = q — p, rename p' as p again, and write T\ as, 



Ti = -if d n xd n yh+ v (x)h + ap (y) j ^e^^K^(p), 



(3.17) 



{P) 4 7 (2tt)« (q 2 - ie) [(p - q) 2 - ie] 

x (rf^ _ p)pqT _ ^ n yP^) ( V Xa ^(q - p) x q c - i{n)p 2 P a ?) . 

The momentum integrals can be computed in the standard way, see appendix D, and 
expanding around n = 4 we get after a rather long calculation, 



■ (n-4) 
15 

+0(n-4) 2 



^8 ( P^ v P a/3 — 3 P M/3 ) +5 P^ v P a/3 ^ 



where 



A(p) 



i / 1 1 , 
r h - In 

8tt 2 U - 4 2 



p — 2e 



+ 0(n-4) , 



(3.18) 



(3.19) 



and the parameter /i is a fixed parameter involving Euler's constant 7. Now using the 
expressions of the Riemann components in terms of the projectors P^ v of appendix B.2 we 
can write, 



a 



1 



/ d 4 x (SR^^R^ix) - R +2 (x)) + d A xR + \x) 



4 In - 4 

- J d 4 xd A y [?>R+ vap {x)R + ^{y) - R + (x)R + (y)\ K,(x - y; /i ) 
+0(ra-4) }, (3.20) 



where 



(p 2 — ie) 



(3.21) 



and a = (28807T 2 ) 1 . Note that the divegent terms with a pole at n — 4, are local and 
quadratic in the curvature. They may be compensated by counterterms in the gravitational 
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part of the action coming from (3.6). Recall that the curvature terms here depend on 
the metric g^ v rather than the physical metric g^ u . 

Let us compute now the fourth term of F$, i.e. T 2 = ^Tr(V^ ] G° + _V a) G°_ + ), which 
depends on the propagators G°^__ = — A + and G°_ + = A~ , represented in (2.22) with 
m — 0. We can write 

T 2 = - l -J d n xd n x'vl 1 \x)A + (x,x')V^\x')A-(x',x) 

= (-«*>) / dVx'j ^^A^WWip ) 



X 



(2vr) n (2tt) 

(d,K(x)) d u + h^(x)d,d u + Z(n)I$\x) 



iq-(x-x') 



ip-(x'-x) 



'(afJ&Wj) d' p + h a - P (x')d'Jp + C(n)R { -\x') 
changing the integration variable from p to q — p as in the previous case we can write, 

Jp-(x-y) 



(3.22) 



T 2 = -8in 2 J d n xd n yh; u (x)h^(y) J 



(27r) n 



IT^ip) = \ J $-6((p - qT)0(q O - P°)S(q 2 )9(-q°) 



(3.23) 



(2 

x (rj pt " TU (q - p) p q T - ^(n)p 2 P >lu ) (rj Xa ^(q - p) x q a - i{n)p 2 P aP ) . 

After performing the phase space integrals and expanding around n — 4, see appendix D, 
we obtain 



where 



J»(p) 



1 



^-^(-P 2 )^(-P°) + 0(n - 4) 

47T 



(3.24) 



(3.25) 



which has no poles at n = 4, therefore the term T 2 requires no counterterms in the action 
to be renormalized. Using the expressions of the Riemann components in terms of P^ (see 
appendix B.2) we get, 

T 2 = || d*xd*y [3R+ ua p(x)R-^(y) - R + (x)R-(y)} K 2 (x - y) + 0(n - 4), (3.26) 

where 



K 2 (x -y) = -\j ^e^-y\2ni)e(- P y(-p ). 



(3.27) 



Here again the Riemann components refer to the metric g pu . We now haver^)=T 1 + T 2 
which must be renormalized by adding the gravitational action up to the second order 
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in hfj, u . The explicit expansion of S g , in (3.6), up to this order in terms of the curvature 
components of the metric g^ u is: 

S g [g^u] = Sg ^ + Sjp + Sg + . . . 



1 



167rGjv 

a 



<fx{-g{x)) 1/2 e^[R{x 
d 4 x[3R^(x)R^(x) - R 2 (x)} 

P 



4(ra - 4) 

+a J d 4 x[R^ a/3 (x)R^(x) - R^(x)R^(x)]\n(^) 

+a J d^xi-gix^ftR^u.^ + RU g oj - A{n gU j)uj, u uj' u - 3(D^) 2 - 2{u. v w 
+0(n-A), (3.28) 

where we have dropped the + sign on the fields for simplicity. Finally adding these terms 
to T m and including only the terms which contribute to the variation of h£ u we get the 
renormalized effective action, 



,(2) 

(«) 



<*,ht]= d'x{-~9 + {x)) l/2 



R+i 



x 



16ttG 



N 



^R + (x)R+(x) 



+2a j d A x(-g + (x)) 1 ' 2 \G + ^ {x)u ;tt u ;v + a g u{u ;v u' v ) + -(^ ; f 
+a ( dH-9 + {x)) l/2 \(RU (x)R + ^(x) - R: u {x)R + ^{x))] u{x) 



+ 



AA(YW) 1/2 (-ff + (!/)) 1/2 

x 3i?+ af3 (x)R + ^(y) - R + (x)R + (y)} K x {x - y; Ji) 
+| J d'xd'y{-g + {x)f/ 2 {-g-{y)f/ 2 



x 



?>R+ va p{x)R-^{y) - R + {x)R-{y)\ K 2 (x - y) 



+0(h%), 



(3.29) 



where fi = /i/io and we have substituted 1 by the volume densities i n au the integrals 
involving quadratic curvature terms of the metric g^ u in order to facilitate the identification 
of the exact variational formulae of the appendix E needed in the computations of the next 
section. 

If one is interested in the production of particles, the in-in effective action is not suitable 
because the probability of particle creation is related to the transition amplitude from the 
in to the out vacua, and this amplitude is directely related to the in-out effective action 
f24| . Let us write the vacuum persistence amplitude in the presence of an external source, 
i.e. (0,out\in,0)j = exp(iW[J}). The probability of pair creation is proportional to the 
imaginary part of W 

P = 2ImW, 
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but if we take J = and consider the quantum fields propagating in the gravitational 
background, then W is just the in-out effective action: r^ ). The calculations leading to 
such an action are similar to those for the in-in case, although they are simpler because we 
do not need to introduce two fields. The main work has already been done: consider just 
a single field in (3.8) and ignore the term Ti- The renormalized action is obtained again 
by adding the gravitational action (3.28), the final result can be read directly from (3.29): 
ignore the plus indices and the term involving K 2 (x — y). The non local term now includes 
Ki(x — y; p) only, this term is complex and it is responsible for the particle creation effect. 
It turns out that the pair creation probability is given, however, by a local term as is well 
known [24]] In fact, the imaginary part of the kernel Ki(x — y; p) is, 



(3.30) 



From the expression for T/ io \ and (3.30), after performing the x and y integrations which 
lead to the Fourier transform of the curvature tensor, R liva p{p)i we get 



an 

P 2/mr (io) = — 



3R,v a p(p)R^ a(3 (- P ) - R{p)R{-p) 



-P 



(3.31) 



Finally, using that the Gauss-Bonnet topological invariant is zero, the above relation can 
be written in terms of the Fourier transform of the Weyl tensor of the physical metric 

as, 

P =9^/(0il^< P »| 2e (- p2 » 9 '- p0 )- < M2 > 
in agreement with the expressions computed by other means [IT1 . 



4 Semiclassical equations 



In this section we obtain the quantum mechanically corrected Einstein's equations due to 
the presence of a massless conformal scalar quantum field. The semiclassical equations 
for the metric perturbation can be found by functional differentiation of the in-in effective 
action (3.29) with respect to h^ v (x) and then restricting \i^ v 



x) 



h^(x) 



as 



$h+ u (x) 



(4.1) 



From equation (4.1) it is easy to derive the equations of motion to first order; we use 
that for an arbitrary functional 



-g8g, 



(4.2) 
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to find the variation of the first two terms of (3.29). Notice that we assume that u(x) is 
a scalar function independent of the metric (in general we will assume that u(x) depends 
on the spacetime point only, in particular, in the flat FRW case it will be a function of the 
cosmological time t only). 

Using the expressions listed in the appendix E one can show that the semiclassical 
equations can be written as 



16nG 



N 



-aR%C^ + ^ [-4(C ( ^), Q , + / #yAft(y)H(x - y; fi) 



+ O(/^)=0, (4.3) 



where G p,v {x) is the Einstein's tensor, C pau/3 (x) the Weyl's tensor, B^ v (x) and A^ u (x) are 
the exact spacetime tensors given by the variation of / d 4 xR 2 (x) and / d 4 xC flotll pC ,w "'P , 
respectively (see appendix E.2 and use that C^pC^ = R^ a upR^ au(3 - 2R llv R iV + \R 2 ), 
with respect to an arbitrary metric g^(x) 

B pv (x) = ^g^R 2 - 2RR" V + 2R pv - 2g^D g R, 

A^(x) = \g pv C a p p(J C a ^ - R^ p R\p p + ARTRJ 

--KBT - 2U g R» v + -RW + -g^n g R, 
3 3 3 



and H^ u (x) is the spacetime tensor defined by 



2 

3' 



H""(,v) = -/?''" /?„" - -JUV U ' f -g^R a pR a P - - A g pv R 2 



1 

V 



In equations (4.3) we are only interested in the expressions of these tensors up to first order 
in the perturbations, as indicated by the bracketed subindices and we recall that an over 
tilde on the tensors refer to the physical metric. The non-local part H(x — y; fi) is the sum 
of the integrals (3.21) and (3.27), 



H(x-y;fi) = K 1 (x-y;fj l ) + K 2 (x-y) 
1 



f d4p e ip - {x ~ y) /in 


(p 2 — ie) 


/ (2tt) 4 \ 


fi 2 



(2ni)6(-p 2 )6(-p° 



(4.4) 



which can be simplified to 



K(x-y;n) = -- J 



d A p 

W) 4 



fi 2 



+ m6{—p 2 )sign{—p Q ) 



(4.5) 



by using ln(±i) = ±i| and lim e ^ m (e + ix) = In \x\ + sign(x)i^. Note that this equation 
is real, in spite of appearences, because the imaginary part of the integrand is an odd term 
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with respect to the integration variable p^. One can also notice that H(x — y; p) differs 
from that defined by Horowitz p8f by a factor (l/Air). 

Equations (4.3) are dynamical equations with higher order derivative terms. When the 
background is flat they reduce to the field equations studied by Horowitz |28| and Jordan 
p0| . As we have noted in the previous section one can add a classical stress matter source 
to these field equations. 

To compare the functional method used in this paper with other techniques and, in 
particular, to compare equations (3.29) and (4.3) with previous results one can give, for 
exemple the energy-momentum tensor of the quantum field and the expression of the 
semiclassical equations in two dimensions. 



4.1 Stress tensor to first order 



From (4.3) one can read the zero and first order vacuum expectation value of the energy- 
momentum tensor of the scalar field, 



TO 



a 



a 



n (o) g^(0) 



1 X 



+3e 



-6lu 



-A{C^u), aP + 



d 4 y^i){y)B.{x - y; fj) 



(4.6) 



(4.7) 



The stress tensor to first order in h^, (Tff(), is in agreement with that obtained by Horowitz 



and Wald and Starobinsky |2"7| . On the other hand, the zeroth order tensor (T/^), 
which gives the exact stress tensor for a conformal scalar field in a conformally flat space- 
time, agrees with that found by other techniques [42] [12]. Note also that we recover the 
trace anomaly result to this order in h^, 



a 



D § R + [R^R 



■fXU 



-R 2 
3 



+ o(h 



4.2 Two dimensional gravity 

As an exercise we will derive here the semiclassical corrections to the stress tensor of the 
scalar field in two spacetime dimensions. At the classical level there is no dynamics for 
the gravitational field in two dimensions because the Einstein tensor is identically zero, 
but the semiclassical corrections due to the presence of a quantum field lead to non trivial 
effects. Since we have worked in n dimensions in order to use dimensional regularization, 
we can use most of the work done in the previous section to analize the two dimensional 
case. From equations (3.17) and appendix D it is not difficult to show that, 

K% D l{q) = ^r 1 ^ - 2)q% u P a , + 0(n - 2) 2 (4.9) 
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where 



ir(q 2 — ie) 



1 7 

+ i + 0(n-2) 



7i-2 2 



(4.10) 



Then, the in-out contribution to the effective action Ti becomes, 

Tt = ~ J d 2 xd 2 yR(x)A (2D \x - y)R(y), (4.11) 

where Ap D ^ (x — y) is the Feynman propagator in two dimensions, 

Ar(l _ s) = _/**^. (4. 12 ) 
J {2iry q A — ie 

From equation (3.23) it is straighforward to see that there is no contribution to the effective 
action due to the typical in- in term T?,. Quantum corrections appear to be non-local and 
quadratic in the scalar curvature, but contrary to the four dimensional case there are 
no divergent terms in the regularization process and counteraction terms are not needed. 
Finally one can express the effective action to one loop as, 



- ±.J d 2 xd 2 y{-g{x)) l / 2 {-g{ y )Y' 2 R{x)A {2D \x-y)R{y) 

+ 0(h%). (4.13) 

The stress tensor of the quantum field can be obtained from the second term of the above 
equation by differentiating with respect to the metric. To order zero in the perturbation 
the energy-momentum tensor vanishes, but to first order there is a non-local contribution, 

<^(*)> (2D) ~ / cfyfrT V - R^]*$ D \x -y) + 0{h%), (4.14) 

and the trace is local {T^ )ti (x)) {2D) = R/(24tt) +0(/i^) as expected [0. Note that in two 
dimensions there is no difference between the in-out and the in-in effective actions, thus 
the use of the usual in-out effective action to derive the semiclassical equations is justified. 



5 Backreaction on the field of a static cosmic string 

In this section, as a simple example, we discuss the backreaction due to one loop quantum 
effects on the gravitational field of a static cosmic string. First, we must compute the 
vacuum expectation value of the stress tensor of a conformally coupled massless scalar 
quantum field outside the core of a straight and static cosmic string. 

In the weak field approximation the metric of a cosmic string can be seen as a small 
perturbation about flat space, 

9fiu(x) = r]^ + h^x). (5.1) 
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Let Tg v be the stress tensor of a cosmic string (or of any other classical source), the 
semiclassical equations up to first order in the perturbation, can be written us, 



G^ix) = 8nG 



N 



Tr(x) + (TlT ) (x))+0(hl)\. (5.2) 



There is no zero order correction, (T^), to the classical Einstein's equations because 
the background is flat; see (4.6). The vacuum expectation value (Tj^) is obtained from 
equation (4.7) when the conformal function is u = 0. We have 

(Tfi(x)) = ~B§{x) +3aJ d*yR(x - y; ji)Aft(y). (5.3) 

where 

Bft(x) = 2r ] ^G a a /-2G a a ^ J 

Aft(x) = -2G""V - \G a a ^ + \rTG a a /. (5.4) 

Note that the quantum correction term (T/Jq) depends on the Einstein tensor G^, thus 
one may use Einstein's equations to the classical order, which is already first order in h^ Vl 
to substitute G^ u by 8ttGnT^ u . This simplifies considerably the problem since the explicit 
gravitational field of the string (or the classical source) is not required to compute (T^). 

The stress tensor of a static cosmic string along the z-axis can be written in the thin 
line approximation as, 

V(x,2/) = -fi5(x)6(y)diag(l,0Al), (5-5) 

where \x is the mass per unit lenght of the string. Since G^fi ~ 10~ 6 for GUT strings we 
can assume that Gat/i <C 1 and the linear approximation (5.1) is justified. In fact, ignoring 
quantum effects the stress tensor (5.5) leads in the linear approximation to the conical 
metric with a deficit angle of 87rGV/i f4"5|| 

Equations (5.4) can be expressed in terms of the string stress tensor and its trace 
T = T C = -2/i5(x)<%), as, 

B [1)lM "(x) = 16ttGW(W -T/), 

A (1) S(x) = 8nG N (-2T„»^ -~T/ + ^T a «). (5.6) 

In this case, B^^ix) in (5.3) gives no contribution outside the core of the string because 
it is proportional to partial derivatives of a delta function with support on the z-axis. The 
tensor components An^ l/ (x), in the Minkowskian coordinates (t, x, y, z) in which equation 
(5.1) is given, are 



v 8ttGn 



/-T xx -T yy \ 

2Tyy 2T X y 

2Ty X 2T XX 

V -T XX ~Tyy 
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The non-local term H(x — y; //) in equation (5.3) can be expressed as a delta function; 
following Jordan |30| we can write 



H(x - y; n) = -—6'[(x - y) 2 } 

Z7T 



(5.7) 



Introducing these expressions into (5.3), we may perform the space integrals; this is easy 
due to the presence of the delta functions. In fact, as an example, let us compute the xy 
component of the stress tensor, (T x y ) = (T y x ), to first order in Gjv/x. The only contribution 
to this component comes from A^ x y = —16TrG^T )Xy /3 = 327iGN^S'(x)S'(y)/3, 



(T/) = 3a / dVH(x - x'; Ji)A {l)x y {x') 
3a 



~ J d 4 x'5'[(x - x') 2 ]A (1)x 



y <x') 



This expression can be written as, 

3a ,. d 



where A(u x v (x ,fi + x M ) is now 



A (1)x y (x'»+xn = 



(5.1 



(5.9) 



(5.10) 



Following straightforward steps, 

d 8 



(T x y ) 



16aGjvA* 



dx dy 



r 99 
lbaGjv/z— — 

ox ay 



lim ^~ 

d 



lim — / dz' 



32aG^ii 



d d 
dx dy 
Axy 



a^o- dX 

oo 



d A x'5[x' 2 - X]5{x' + x)5{y' + y) 
1 



^/z' 2 + x' 2 + y' 2 — X 



dz' 



o {z 12 + x 2 + y 2 ) 3 / 2 



(x 2 + y 



2\3 



(5.11) 



The remaining non null components of the stress tensor in Minkowskian coordinates to the 
same order GnU are, 



(T y y ) = 32aG N fx 
(T/) = 32aG N fi 



y 2 - 3x 2 
(x 2 + y 2 ) 3 
x 2 - 3y 2 



x 2 + y 2 ) 3 
T t l ) = (T z z ) = 32aGV/i 



(x 2 + y 2 ) 2 



(5.12) 
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Because of the cylindrical symmetry of the problem it is better to express this tensor 
components in polar coordinates as, 

<V> = ^^^(1,1,-3,1), (5.13) 

which is independent of the arbitrary renormalization scale ft, (as expected because it 
would contribute with delta functions with support on the core of the string which we do 
not consider). This tensor coincides with the first order development in G^n of previous 
exact results obtained by other techniques f37|. Note that in the exact case, i.e. when the 



classical solution is found explicitly, the one loop quantum stress tensor is simply (5.13) 
where one changes 32GnH by 2[(1 — 4Gnh)~ 4 ~ !]• 

The backreaction equations (5.2) have been solved by Hiscock [R!J who found that 
the linear corrections to the metric outside the string are such that the spacetime is no 
longer flat space with a deficit angle: the two surface perpendicular to the string is an 
hyperboloid which asymptotically approaches the conical surface at large distances (the one 
loop quantum corrections to h^ v {x) are of the form Gn^Ti/v 2 ). Note that the semiclassical 
equations here have no higher order derivatives because we have treated the quantum terms 
as a perturbative correction (as Hiscock does) in line with Simon's arguments [32]. Work 
on the backreaction on dynamic cosmic strings |||J is in progress. 
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7 Appendix 



A 



. Useful relations 



A.l. Bianchi identities 



j^fiaul3;a _|_ jyiafics,v _|_ R tiacru;f3 _ g (A 1) 

2R a P. a = R'P. (A.3) 
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A. 2. Commutation of covariant derivatives 



k=l 



R^\ = g ua R^. pa -R^R ap + R^R a v 



R mVP -0a = R» aUl3 ]a( 3 = n g R^-R^ a = 

= DgRT - -R^ v + R» aul3 R a(3 - R» a R a 



A. 3. 2D curvature tensors 



R 



Gfxv = 0. 



t^iiavfit^ — It — Iti^ti . 



B . Expansions around flat space 
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.1. Curvature tensors 



g^u(x) =r] IMU + h^(x), (-,+,•••,+). 



g ^( x ) = jf - h^(x) + h\{x)h au {x) + 0{hl v ). 



-g{*)) 1/2 = 1 + l h + l h2 - \ h ^ v + o(h%). 



pa a pa a pa . pa pM pa pM 

-rt /3-yS — 1 l gs — O5I p 1 -+- 1 M7 1 ^ — 1 m5 1 ^ 



f3~ ( 8 — /3«5 

1 I 



S\5,f3 = h\s,{3 + hp X ,8 — h/3S,\- 



(B.l) 
(B.2) 

(B.3) 
(B.4) 

(B.5) 
(B.6) 
(B.7) 



R 



05 



flUJ 



(B.8) 



R — h a ^ :a p — h :a a 

3 1 

-\h^ a K^ - ibrh^s + h^h pu + rv; + o(h 



R, Pl sR aM& = Ks^h" 5 ^ - 2h aS ^h^ + Ks^h?™* + 0(h 



3 ) 



(B.9) 
(B.10) 



RpsR 135 — ^ 



+ ^,a° V/ + 2^,a" V + ^ V + O(^) 1 . (B.ll) 
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R 2 = h a ^ aP hr,^ - 2h a P, aP h,f + h a a h/ + 0(h%,). (B.12) 



B.2. Curvature tensors in terms of the projector P t 



(B.13) 



G = J d n x (R^pR^ - AR^R^ + R 2 )=0 + 0(h%,). (B.14) 
J d n xR 2 (x) = J d n xd n yh^(x)h a p(y) J ^-e^^P^P^q 4 + 0(h%). (B.15) 



/ d n xR^(x)R^(x) = J d n xd n yh, u {x)h aP {y) J *^ e *<*-v) P tf P ™ q 4 + 0{h ^y 



(B.16) 



J d n xd n yR(x)R(y)K 1 (x - y; (i) = -- J d n xd n yh liu (x)h a p(y) 

P^P a Pq* + 0(h%,). (B.17) 



x , '■ In 

{2n) n 



q — le 



j dVyR^WR^iy^x-y-ii) = ~ J d n xd n yh, v (x) h af) (y) 



x f In 

(2vr) n 



q — it 



pnPpva q 4 _|_ Q(h 3 ). (B.18) 



Jd n xd n yR(x)R(y)K 2 (x-y) = - l - J ' d n xd n yh^{x)h a p{y) 



x / ^e^-y\2m)e(-q 2 )9(-q°)P^P^q 4 + 0(hl l/ ). (B.19) 



(27T) 



Jd n xd n yR^(x)R^(y)K 2 (x-y) = ~ j ' d n xd n yh, u (x)K p (y) 



I ^** <X ^™W-<?M-f)^^Q A + 0{h%). (B.20) 



(27T) 
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C Curvature tensors in conformally transformed n diamensional spaces 

9 ltv {x) = e 7u g^{x) (C.l) 

-1g a [ngi3v]U-pUJ' p ] . (C.2) 

R/jot = R^ a - in- 2)u ]lia - g m O g uj + (n - 2)u ;ii u ;a - (n - 2)g lia u).^ p . (C.3) 

R = e ~ 2uJ [R - 2(n - l)U g u - (n - l)(n - 2)uj, p uj'' p ] . (C.4) 

R, va pR pvaP = e" 4w [R^R pua/3 - 8R^(u.^ - u,^) - ARw,^ 

+4(n - 2)uj.p V u)'^ v - S(n - 2)u' pv uj. p uj ]V + 4(D^) 2 

+8(n - 2)uj.^uj'' p U g uj + 2(n - l)(ra - 2){uj, p uj^) 2 ] . (C.5) 



R^BT = e' 4uJ \R^R^ - 2{n - 2)R^(u ;tlu - u,^u, v ) - 2RU g u 

-2(n - 2)Ruj. ix uj^ + (n - 2) 2 uj.^ v + (3ra - 4)(D g cu) 2 
-2(n - 2)V'"V /i u; ;i , + 2(n - 2)(2n - 3)cj ; ^ ; ^D 9 c«; 
+(n-l)(n-2) 2 (o; ; ^) 2 ]. (C.6) 



^ 2 = e~ 4cJ [i? 2 - 4(ra - l)RU g uj - 2(n - \){n - 2)Rw.^ + 4(n - l) 2 (n^) 2 

+4(n - l) 2 (n - 2)uu.^n gU o + (n - l) 2 (n - 2) 2 {u,^) 2 } . (C.7) 



4 ~ ~ v 2 

=~2) RvvBr + (n-l)( 5 
= e^C^C^ (C.8) 



p p u^vafB _ p p^f j p2 

°^ Q/3 ° - ( n -2) ^ (n-l)(ra-2) 



D . Momentum integrals and dimensional regularization 
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d n p 



f(p,q)- 



(2tt)' 

d n p £l N I(q) 



1 fll/ 



h = J j2^fMP* = -2 



—i p 

(n-1) 



(27T) 

m 



-fip, q)P»PvPa 



(n-1) 



/d n p 
J^y^f{P,Q)P^PuPaPf3 

I(q) f q 2 

= QnVvQaVP — 



(n-1) 



(n 2 - 1) 



[P/j.uPa/3 + PuflPva + PvfiP\xo\ 



h( P ) = J 



d n q 



(2tt)™ (q 2 -ie)[(p-q) 2 -it] 



— i(p 2 — ^e) 2 



; r( 2 -f) [r(f-i) 



(47r)fr(n-2) 



D) 



(P) = 



8tt 2 



_1_ + I ln( p 2 - i e ) + I ( 7 - 2 - In 4tt) + C>(n - 4) 



li 2D \p) 



7r(p 2 — ie) 



+ l -l + 0{n-2) 



n-2 2 



hip) 



d n q 



5(q 2 )e(-q°)5[(p - q) 2 ]6(q° - p°) 
1 



(2tt / 

9(-p°)9(-p 2 )(p 2 ) n ~ 3 



{n-3)2 n - 1 (2n) 3 \p\ 



1 



(P° + Ip1) n_3 (P°-\p\) n ~ 3 
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r (4D) 
-to 



(P) 



4 2D) (p) 



g(-p°)g(-p 2 ) 

4(2vr) 3 ' 

e(-p°)e(- P 2 ) 

(2tt) 3 (p 2 - ie) 



E. Variational calculus 
E.l. Variational equations 



i 

ST ps = ~g aX [$gx8;p + $gp\-,s - Sg P s;x\ 



I 

5R = -R ps 8g p5 + /V A [Sgx8;p a - Sg mXa ] . 
5(n g u) = -uj'^Sg^ - -u'' x g' tv [Sg^ + Sg^ X;u - Sg^x] . 



E.2. Functional differentiation 



5 J d 4 x^g~R 2 = J d 4 xy/^g^-g^R 2 - 2RR^ + 2R^ - 2g^n g R^ 5g 
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5 J d 4 Xyf=g~R al3 R a p = J d 4 x^^g^R a PR ap - 2R^R a v + 2R^ V a - U g R^ 

~g^a g R^Sg^. (E.10) 

5 J d 4 x^g-R a ^R aPpa = J d 4 x^^g^R a ^R af3pa -2R^R\ (5p -AU g R^ 

+2R'" U - AR» avp R af} + AR» a R a v } Sg^. (E.ll) 



5 J d 4 x^R a(S to. a u ](3 = J d 4 x^=g~ | g^R ap u, a u, p - 2R^ a u' v) u- a + {U g u)u^ v 



Q 

2 y 



U g uf + u'^u.^ + [co^ a - <ru g (con) co ;a \ 6g^. (E.12) 



5 / d 4 x^~gR (D g u) = / d 4 iv^ { -R^ (n g cu) + R^u^ - -g^R a u ]a + (□, 



,UJ) 



-g^U g (n g u) 



(E.13) 



5 J d 4 x^=g~R (uj. a uj' a ) = J d 4 x^ 



^g- l^-g^R (cu ]a u-> a ) - (u^) - R (uV) 
+2R^cu. >a cu.p + 2 



+2 [u^ a - <T D> ;a )] u, a } 5g, v . (E.14) 



5 J d 4 x^=g~U g u (co. a ^ a )=J d*xy/=g feu^u^ - g^u^u.j 



CO 



a - {D gU j) u V J Sg^. 

(E.15) 



5 J d 4 x^g-(U g u f = J d 4 xV=g 



(E.16) 



5 J d 4 x^Tg-{uj, a ^ f = J d 4 x^^-g^{u, a w a f - 2uV (u;. a u;' a )} Sg^. (E.17) 
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5 J d A x % f=g~R t * v R IJkV oj{x) 



: J d 4 x^g { \^-g^R^R a p - 2R^R aP - R^ 

+ r " 



+2K*bu>% - R»» (O g Lo) - g^R aP u. aP } 5g, u . (E.18) 



5 J d 4 x^gR a ^R af3p Mx)= I d A x^-g\^g^R a ^R aPpa -2R^R 



-4R» auf3 R a p + ±RTR a v - 4n g RV + 2R'" U ] uj{x) 
-8 [R^' a - R a ^] cu ]a - AR^u ;afi } 5g pu . (E.19) 

= J d 4 x^-g+(x) {2 J d%^-g+(y)R + (y) 

-G + ^(x) + VfoVfo - g + ^(x)D {x) } K ± (x -y-p)\ Sg^x). (E.20) 



5 J d 4 xd^-g+(x)^/-g+(y)R +a ^(x)R^ pa (y)K 1 (x - y; fx) = 

= J (Pxy/-g+(x) [J d^-g+(y) [g + ^(x)R^ pa (x)R +a ^(y) 

-4^(y)VWvf] Mx -y-Jx)) 5g+,(x). 



(E.21) 



5 J d i xd i y^-g+(x)^-g-(y)R + (x)R-(y)K 2 (x-y) = 

= J d*xy/-g+(x)y d i y^-g-(y)R~(y) 

-G + ^{x) + VfoVfo - g + ^(x)D {x) ] K 2 (x - y) \ Sg^x). (E.22) 



5 J d 4 xd%^-g+(x)^-g-(y)R^ pa (x)R- a ^(y)K 2 (x-y) = 

= J d*xy/-g+(x){J d*yy/-g-(y)] \^9 +lxv {x)R^ pa {x)R-^ \y) 

-2R-^{y)V^Vf 1 K 2 (x - y) ) 5g+,(x). (E.23) 
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